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Abstract. We introduce the notion of the descent set polynomial as an alternative way 
of encoding the sizes of descent classes of permutations. Descent set polynomials exhibit 
interesting factorization patterns. We explore the question of when particular cyclotomic 
factors divide these polynomials. As an instance we deduce that the proportion of odd 
entries in the descent set statistics in the symmetric group 6 n only depends on the number 
on l's in the binary expansion of n. We observe similar properties for the signed descent set 
statistics. 



1. Introduction 

The study of the behavior of the descent sets of permutations in the symmetric group 6„ 
on n elements usually involves such questions as maximizing the descent set or determining 
inequalities which hold among the entries [31 [HI El EH 031 EH EES] ■ The usual way to encode the 
descent statistic information is via the Eulerian polynomial A n (t) = ^2 S (3 n (S) ■ where S 

runs over all subsets of [n — 1] = {1, . . . , n — 1}, and (3 n {S) denotes the number of permutations 
of size n with descent set S. We instead introduce the descent set polynomial where the statistic 
of interest appears in the exponent of the variable t, rather than as a coefficient. That is, the 
nth descent set polynomial is defined by 

o„(*) = E* A,(fl) . 

s 

where S ranges over all subsets of [n — 1]. 

The degree of the descent set polynomial is given by the nth Euler number, which grows 
faster than an exponential. Despite this, these polynomials appear to have curious factorization 
properties, in particular, having factors which are cyclotomic polynomials; see Table [5] This 
paper explains the occurrence of certain cyclotomic factors. We have displayed these in boldface 
in the tables. Both combinatorial and number-theoretic properties (for example, the number of 
l's in binary expansion of n and the prime factorization of n) are involved in our investigations. 

The divisibility by cyclotomic factors is related to the remainders of sizes of descent classes 
modulo certain integers. As a simplest example, Q n (t) is divisible by the second cyclotomic 
polynomial <j> 2 if and only if the number of even descent set classes is equal to the number of 
odd descent set classes. In other words, the proportion of even and odd entries in the descent 
set statistics is the same (in the notation below, p(n) = 1/2) if and only if —1 is a root of the 
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Table 1. The proportion p(n) for at most five l's in the binary expansion of n. 



descent set polynomial. Somewhat surprisingly, whether or not n has this property depends 
only on the number of l's in the binary expansion of n. 

The paper proceeds as follows. In Section [2] we look at the proportion of odd entries in the 
descent set statistics. In Section [3] we discuss this result from the viewpoint of quasisymmetric 
functions related to posets. We consider similar properties for the signed descent set statistics 
in Section |H The natural setting for this question is to look at flag vectors of zonotopes. In 
Section [5] wc explore patterns of descent statistics modulo 2p for prime p. Here we introduce 
the descent set polynomial and consider divisibility by cyclotomic polynomials. In Section [7] 
we explore when the descent set polynomial is divisible by the quadratic factors $2, $4 and 
&2 P - I n Section [5] we introduce type B quasisymmetric functions and the signed descent set 
polynomials. We use the former to describe divisibility patterns of the latter. Finally, in the 
concluding remarks we make a number of observations on the data presented in Tables [5] and [3J 



2. The proportion of odd entries 



For 7r — 7Ti ■ ■ • 7r„ a permutation in 6„, recall that the descent set of 7r is the subset of [n — 1] 
given by {i : iti > 7T; + i}. For a subset S of [n — 1] the number of permutations in 6 n with 
descent set S is denoted by (3 n (S). 

Let p(n) denote the proportion of odd entries in the descent statistics in the symmetric 
group 6„, that is, 

( , _ \{SC [ n -l] : = 1 mod 2}| 

P\ n > ~ 2™- 1 
For instance, p(3) = 1/2 since in the data 3 (0) = /3 3 ({1,2}) = 1 and /3 3 ({1}) = /3 3 ({2}) = 2 
exactly half of the entries are odd. 

The first few values of the proportion p(n) are shown in Table [T] In this section we prove 
the following result. 

Theorem 2.1. The proportion of odd entries in the descent set statistics p{n) only depend on 
the number of 1 's in the binary expansion of the integer n. 

Recall a composition of n is a list 7 = (71,72, ■ ■ ■ ,7m) of positive integers such that 71 + 
72 + • ■ ■ + 7 m = n. The multinomial coefficient is defined by 

/ n\ n\ 

VfJ 7i ! • 72! • • -7m!' 
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Define a bijection D between subsets of the set [n — 1] and compositions of n by sending the 
set {si < S2 < ■ ■ ■ < s m -i} to the composition (si, S2 — si, S3 — S2, . . . , n — s TO _i). Let a n (S) 
denote the multinomial coefficient ( D %y) ■ The following is a classic result due to MacMahon: 

Lemma 2.2. Let S be a subset of [n— 1]. Tften £/ie number of permutations in & n with descent 
set contained in S is a n (S), and we have 

Pn{S)= £(-l) |S - T| -anCr). 
TCS 



We need Kummer's theorem for the multinomial coefficient version. 

Theorem 2.3. For a prime p and a composition 7 = (71,72, ■ ■ • ,7m) of n, the largest power 
d such that p d divides the multinomial coefficient (") is equal to the number of carries when 
adding 71 + 72 + • • • + 7m in base p. 

As a corollary we can determine whether a multinomial coefficient is even or odd. This 
corollary also follows from Lucas' congruence for binomial coefficients. 

Corollary 2.4. For a composition 7 = (71,72, ■ ■ ■ ,7m) of n, the multinomial coefficient (™j is 
odd if and only if there are no carries when adding 71 + 72 + ■ • ■ + 7m in base 2, that is, for all 
i =/= j , the binary expansions of 7, and jj have no powers of 2 in common. 

Let the binary expansion of n be n = 2 n + 2 J2 + • • • + 2^ k , where ji > J2 > • • • > jk ■ Call an 
element of [n — 1] essential if it can be expressed as X^es f° r some nonempty proper subset 
B of [k]; otherwise, call this element nonessential. 

Lemma 2.5. If S C [n — 1] contains a nonessential element Si, then a n (S) is even, that is, 
a n (S) = mod 2. 

Proof. Let 7 = (71, 72, • ■ ■ , 7m) be the associated composition D(S). Notice in the addition 
(71 + • • • + 7i) + (li+i + • ■ • + 7m) = Si + (n — Si) = n there is a carry in base 2. Hence it follows 
from Corollary 12.41 that a n (S) is even. □ 

Lemma 2.6. Let S be a subset of [n — 1], and suppose that i 6 [n — 1] — S is an nonessential 
element. Then 

f3 n (S)=f3 n (SU{i}) mod 2. 



Proof: By Lemmas 12.21 and 12.51 we have 

/J(su{i}) = ^(-lj'^-^m + ^Hpi.^truW) 

TCS TCS 

= -/?„(£) +^(-l)l s - T l-a„(TU{ l }) 

TCS 

= /3„(5)mod2. □ 
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Lemma 2.7. Let S = {si < s 2 < • • • < s m _i} fee a subset of [n — 1] consisting of essential 
elements, so that there are nonempty proper subsets B-l, _B 2 , ■ ■ ■ , B m -i of [k] such that s r = 
^2 ieB 2 Ji . Then a n (S) is odd if and only if B\ C _B 2 C • • ■ C _B m _ 1 . 

Proof. Let Bq = and B m = [fc]. If Bi C B 2 C ■ • • C B m _i then 7,, = s r — s r _i = 
Xligs -b 1 ^ ■ Then there is no carry in the addition 71 + • • • + 7, n = n and hence a n (S) is 
odd. On the other hand, if a n (S) is odd then there is no carries in the addition 71 + ■ • • +7, n = n, 
so all the 2-powers that appears in 71, ... , j m must be disjoint. Since s r is given by the partial 
sum 71 + • • -+7r, the 2-powers appearing in s r must be contained among the 2-powers appearing 
in s r +i, that is, B r C B r+ \. □ 

Lemma 2.8. Let n have k \ 's in its binary expansion. Let E = {ei, e2, . . . , e 2 fc_i} be the set 
of essential elements of [n — 1], where the Cj 's are listed in increasing order: e\ < e 2 < • • ■ < 
e 2 h_i- Let S = {si l; s, 2 , . . . , Si m } fee a subset of E and S be the set of indices of S, that is, 
S = {ii, ii, . . • , i m }. Then the parity of /3 n {S) is the same as the parity of f3 2 k_ 1 (S). 

Proof: From Lemma 12.71 it follows that the parity of a n (S) is the same as the parity of 
a 2 k -i(S)- Now the result follows by 

Pn(S)= ^(-l) |s - T| -a„(T)= ^(-l)l^ f l-a 2fc _ 1 (f) = /3 2 ,_ 1 (^)mod2. □ 

TCS fc§ 

We are now ready to prove Theorem 12. II 

Proof of Theorem I2.lt It follows from Lemma [2T61 that the proportion of odd entries among 
Pn(S) for SC. [n — 1] is the same as the proportion of odd entries among (3 n (S) for S C E. 
But by Lemma \2.8\ the proportion of odd entries among (3 n {S) for S C E depends on k, the 
number of l's in the binary expansion of n. □ 



3. QUASISYMMETRIC FUNCTIONS AND POSETS 

In this section we relate the preceding result to the theory of quasisymmetric functions. 

Consider the ring Z[[iui, io 2 , . . .]] of power series with bounded degree. A function / in this 
ring is called quasisymmetric if for any sequence of positive integers 71,72, ... , 7 m we have 

Iw] 1 ■■■w7™] f= \w] 1 ---w^]f 

L H 2fc J ■> V 31 3k 1 ■> 

whenever »!<•••< i m and ]!<■■■< j m , and where [u> 7 ]/ denotes the coefficient of w 1 in /. 
Denote by QSym C Z[[u>i,«; 2 , . . .]] the ring of quasisymmetric functions. 

For a composition 7 = (71,72, . . . ,7m) the monomial quasisymmetric function M 7 is given 

by 
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Definition 3.1. Let P be a graded poset with rank function p. Define the quasisymmetric 
function F(P) of the poset P by 

f(p) = J2m p(c) , 

c 

where the sum ranges over all chains c — {0 = xq < x\ < ■ ■ ■ < x m = 1} in P, p(c) denotes 
the composition (p(xo, X\), p{xi, x%), ■ ■ ■ , p(x m -i, x m )), and p(x,y) denotes the rank difference 
p(x,y) = p(y) - p(x). 

The quasisymmetric function of a poset is multiplicative, that is, for two graded posets P 
and Q the quasisymmetric function of their Cartesian product is given by the product of the 
respective quasisymmetric functions 

F(P x Q) = F(P) ■ F(Q); 

see Proposition 4.4 in [6]. Recall that the Boolean algebra B n is the Cartesian power of the 
chain of two elements B\, that is, B n = B\ and F(B\) = w\ + W2 + • ■ • . Hence we have that 
F(B n ) = F(Bi) n = Mfc =( Wl +w 2 + --- )". 

Let P be a graded poset of rank n. For a subset S = {s\ < S2 < ■ ■ ■ < s m _i} of [n — 1], 
define the flag f -vector entry fs to be the number of chains {0 = xq < x\ < ■ ■ ■ < x m = 1} 
such that p(xi) = Sj for 1 < i < m — 1. The flag h-vector is defined by the invcrtiblc relation 

%=E(- 1 ) |S " T| -^- 

TCS 

Recall the bijection D between subsets of [n — 1] and compositions of n defined in Section [2 
By abuse of notation we will write Ms instead of Mjj^g^ where the degree of the quasisymmetric 
function is understood. Then Lemma 4.2 in [5] states that the quasisymmetric function of a 
poset encodes the flag /-vector, that is, 

F ( p )= 2 fs-Ms. 

SC[n-l] 

The fundamental quasisymmetric function Lt is given by 

L T = M S- 

TCS 

Note that one can write F(P) in terms of fundamental quasisymmetric functions: 

F ( p )= h s-L S - 

SC[n-l] 

If the poset P is a Boolean algebra B n then hs = f3 n (S). This is straightforward to observe 
using the classical i?-labeling of the Boolean algebra (TTl Section 3.13] or by direct enumera- 
tion p~7l Corollary 3.12.2]. The multiplicative property Bi+ m = Bt x B m allows us to compute 
F(B n ) modulo 2: 

Lemma 3.2. Form = 2 J we have F(B m ) = Mi m ) mod 2. Consequently, forn = 2-' 1 +• • ■ + 2 : > k 
with ji > ■ ■ ■ > jk > we have F(B„) = Y\i—i -^f(2«) m °d 2. 
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Proof: It is enough to prove the first statement. Recall the congruence (a + b) 2J = a? 3 + 
b 2 ' mod 2. Now we have 

F{B 2J ) = (wi + w 2 + ■ ■ • ) 2J = wf +wf +■■■ = M (2 i) mod 2. □ 



An ordered partition ir of a set [k] is a list of non-empty pairwise disjoint sets (B±, B 2l ■ ■ ■ , Bj) 
such that their union is [k\. 

Theorem 3.3. For positive integers mi, m 2 , . . . , m^, the product of the quasisymmetric func- 
tions M( mi ) • M( m2 ) • ■ ■ M (mk) is given by 

M {mi yM {m2 yM {mk) =Y, M SM 

where 

S(ir) = L> _1 ^ m t , ^ m„ . . . , m, 
and 7r = (Si, B 2 , ■ ■ ■ , Bj) ranges over all ordered partitions of the set [k]. 

Proof. The theorem can be proved by iterating Lemma 3.3 in [B]. □ 

In view of Lemma 13. 2\ the following is a restatement of Theorem 12.11 in the language of 
quasisymmetric functions: 

Theorem 3.4. The proportion of odd coefficients in the quasisymmetric function f = A/( 2 ji) ■ 
M( 2 J2) • • ■ -W(2Jfc) when expressed in the L-basis only depends on k. 

Proof. There is a natural partially ordered set life on ordered partitions of [k] with the cover 
relation tt -< a whenever a is obtained from 7r by merging two adjacent blocks. Expressing the 
quasisymmetric function / in terms of the L-basis, we get 

/ = J2 M SM 

7ren fc 

= E E (-i) |t - swi -£t, 

where S(ir) is defined as in Theorem 13.31 for jjjj = 2 Ji . Thus the coefficient of Lt is given by 
the sum 

J2 (-l)l r -*«l. 

7ren fc : S(tt)CT 

As in the proof of Theorem 12.11 considering pairs of sets T and T U {i}, where i T is a 
nonessential element, we conclude that the proportion of odd coefficients in the expansion of / 
in the L-basis only depends on sets T consisting solely of essential elements. The coefficients 
corresponding to such T depend only on the poset Hfc, that is, only on k, because of the above 
expression for the coefficient of Lt- □ 

Theorem 12.11 implies that for n = 2- 7 all the entries in the descent set statistics are odd. 
Hence it is interesting to look at this data modulo 4. 
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Theorem 3.5. For n = 2 J > 4 exactly half of the descent set statistics are congruent to 1 
modulo 4, and the other half are congruent to 3 modulo 4. 

Proof. First we claim that F(B n ) = Mr n ) + 2Mr n /2, n /2) mod 4. This identity follows from the 
observation (a + 2b) 2 = a 2 mod 4 and by induction on j. where the induction step is 

F(B 2n ) = (M (n) + 2M ( „ /2i „ /2) ) 2 = M 2 n) = M {2n) + 2M (n>n) mod 4. 

We have the expansion 

M (n) + 2M (n/2>n/2) = M + 2M { „ /2} 

= ^(-i)i s i.L S - 2 . 2 

S n/2SS 

= £ (-i)I*i.l 5 - £ (-1)' S '-L S . 

Hence the descent set statistics modulo 4 are given by /3 n (S) = (— 1)I S ~W 2 }I mod 4. Thus 
for 1 ^ S the values of (i n {S) and ^(S 1 U {1}) have the opposite sign modulo 4, proving the 
result. □ 



4. The signed descent set statistics 



A signed permutation of size n is of the form rr = 7Ti • • • n n where each m belongs to the set 
{±1, . . . , ±n} and 7Ti| ■ ■ ■ |7r„| is a permutation. Let 6^ be the set of signed permutations of 
size n. For ease of notation put tto = 0. The descent set of a signed permutation tt is a subset 
of [n] defined as {i : 7r.;_i > m}. For S C [n] let (3^(3) denote the number of permutations 
in with descent set S. 

An equivalent way to using quasisymmetric functions to encode the flag /-vector data of a 
poset is via the ab-index. Let a and b be two non-commutative variables. For S C [n — 1] let 
us be the monomial u\Ui ■ ■ ■ u„_i where Wj = a if i ^ S* and = b if i G 5. The ab-mdex of 
a poset P of rank n is defined as the sum 

s 

When the poset P is Eulerian then its ab-index can be written in terms of c = a + b and 
d = ab + ba. This more compact form removes all linear redundancies among the flag vector 
entries [2j. The linear relations satisfied by the flag /-vectors of Eulerian posets are known 
as the generalized Dehn-Sommerville relations pQ. Similarly to quasisymmetric functions, the 
ab-index and cd-index also have an underlying coalgebra structure. For more details, see 

The poset associated to signed permutations is the cubical lattice C„, that is, the face lattice 
of an n-dimcnsional cube. Observe that C n has rank n + 1. We have 

SC[n] 

A more general setting for the cd-index of the cube is that of zonotopes. Recall that a 
zonotope is a Minkowski sum of line segments. Associated to every zonotope Z there is a 
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central hyperplane arrangement 7i. Let L be the intersection lattice of the arrangement Ji. 
A result by Billera, Ehrenborg, and Readdy [3] shows how to compute the cd-index of the 
zonotopc from the ab-index of the intersection lattice L. First, introduce the linear map to 
from Z(a, b) to Z(c, d) defined on an ab-monomial as follows. Replace each occurrence of 
ab by 2d and then replace the remaining letters by c. The main result in [3j states that the 
cd-index of the zonotopc Z is given by 

(4.1) $(Z) =w(a •#(£))• 
In particular, for the cubical lattice we have 

(4.2) #(<?„)= w(a.tf(B n )), 

since the associated hyperplane arrangement is the coordinate arrangement and its intersection 
lattice is the Boolean algebra. 

Considering Equation (|4.1[) modulo 2, we observe that *B(Z) = c™ mod 2 and hence we obtain 
the following result. 

Lemma 4.1. All the entries of the flag h-vector of a zonotope are odd. In particular, all the 
signed descent set statistics are odd. 

In order to understand the flag /i-vector modulo 4, we need a few lemmas. 
Lemma 4.2. After expanding the cd-polynomial 

c l ■ d • c"-*" 2 

i=0 

into an ab-polynomial, exactly half of the coefficients are odd. 



Proof. Since d = ab — ba mod 2, it is sufficient to consider the identity 

n-2 

c* • (ab - ba) ■ c™^" 2 = a • (a + b) n ~ 2 • b b ■ (a + b)"~ 2 • a. 

i=0 

This identity holds since the coefficient of an ab-polynomial in the sum is the number of 
occurrences of ab minus the number of occurrences of ba in the monomial. This difference only 
depends on the first and last letter in the monomial and the identity follows. To complete the 
proof, observe that out of 2" ab- monomials of degree n, exactly 2™ _1 appear in the right-hand 
side of the identity. □ 

Lemma 4.3. Let z and w be two homogeneous polynomials in Z(a, b) of degree m and n, 
respectively, each having exactly half of their coefficients odd. Then the two ah -polynomials 

c l ■ z ■ c J and z ■ c™ + c m • w 

also each have exactly half of their coefficients odd. 



Proof. We only prove the second statement of the lemma. Wc omit the proof of the first, as it 
is similar and easier. Let u and v be two ab-monomials of degrees m and n, respectively. The 
coefficient of u ■ v in z ■ c" + c m • w is given by the sum of the coefficients of u in z and of v in w. 
Hence the coefficient of u ■ v is even when the coefficients of u and v are both even (2 m_1 • 2 n_1 
cases) or the coefficients of u and v are both odd (2 m_1 • 2 n_1 cases). □ 



CYCLOTOMIC FACTORS OF THE DESCENT SET POLYNOMIAL 



Combining Lemmas 14.21 and 14.31 we have: 
Proposition 4.4. Let ay, . . . , a„ be integers, not all of which are even. When the cd-polynomial 

n-2 

oh ■ c l ■ d ■ c"" 1 " 2 

i=0 

is expanded into an ah -polynomial, exactly half of the coefficients are odd. 

Theorem 4.5. For a zonotope Z either (i) exactly half of the flag h-vector entries are congruent 
to 1 modulo 4, and the other half are congruent to 3 modulo 4; or (ii) all the flag h-vector entries 
are congruent to 1 modulo 4. 

Proof. Considering the identity (|4.ip . we observe that the only terms in the right-hand side 
with non-zero coefficients modulo 4 are c" and those cd-monomials having exactly one d, that 
is, 

$(Z) = c" + 2 • ( JTj a, • c l • d • c"-'- 2 ^ mod 4. 

If all the aj's are even then the flag h- vector entries are congruent to 1 modulo 4. If at least 
one a.i is odd, then by Proposition 14.41 exactly half of the flag h- vector entries are congruent to 
1 modulo 4 and the other half arc congruent to 3 modulo 4. □ 

Now we consider the cubical lattice, that is, the signed descent set statistics modulo 4. 

Theorem 4.6. For an integer n > 2, exactly half of the signed descent set statistics are 
congruent to 1 modulo 4, and the other half are congruent to 3 modulo 4. 

Proof. Observe that there are 2 n atoms in the cubical lattice. Hence hiy\ = 2" — 1 = 3 mod 4. 
Thus the result follows from Theorem 14.51 □ 



5. Descent set statistics modulo 2p 

For a set S of integers and a non-zero integer q, define the two notions q ■ S and S/q by 

q ■ S = {q ■ s : s <E S}, 
S/q = {s/q : s£S and q | s}. 

Recall that a n (S) (resp., f3 n (S)) denotes the number of permutations in ©„ with descent set 
contained in (resp., equal to) S. 

Proposition 5.1. Let q = p , where p is a prime and t is a non-negative integer. Let n — r ■ q, 
where r is a positive integer. Then the descent set statistics modulo p are given by 

/3„(S) ee (-1)I^[-HI . fos/q) modp, 

where SC. [n — 1] . 
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Proof. Observe that 

M{ 1) =F(B r )= Y, <*r(S)-M s . 

SC[r-l] 

In this quasisymmetric function identity make the substitution Wi i — > wf. We then obtain 
M U = E a r {S)-M q . s 

SfC[r- 1] 

= E E ms) ■ • l t 

SfCfr- 1] rc|»-i] 

_ ,-SCT 

= E E «r(5). (-1)^1. Lt 
rc[n-l] 9-SCT 

= E (-I) |r -«' IP - I]l • f E «r(«) ' (-l) |T/ "- S| ] • LT 
TC[n-l] \SCT/g / 

= E (-i) |r - 9,[r - 1]l -A-(r/g).Lr. 

TC[n-l] 

Since Mf = M( g ) mod p, we have F(B n ) = (Mf) r = MT-. mod p. Now by reading off the 
coefficients of Lt, the result follows. □ 

Corollary 5.2. Let q = p , where p is a prime and t is a non-negative integer. Then 

p q {S) = {-l)^modp. 



Proof. The claim can be deduced from Proposition 15.11 by setting r = 1. A direct argument 
proceeds as follows. Since (a + b) q = a q + b q mod p, we have 



F(B q ) = (w 1 +w 2 + ---) q 
= w\ + w\ + ■ ■ ■ 

= M (i) =E^ 1 )' 5 ' mod p. 



□ 

Corollary 5.3. Let q = p , where p is a prime and t is a non-negative integer. Then 



f3 2q {S) = (-l)l s -{«>l mod p. 



Proof. Follows from Proposition 15 . 1 1 by setting r — 2 and noting that 2 {0) = /^({l}) = 1. □ 



For n having the binary expansion n = 2 J1 + 2 J2 + • • • + 2 Jfc , where j\ > j 2 > ■ ■ ■ > jk > 0, 
recall that an element j € [n — 1] is nonessential if j is not a sum of a subset of {2- J1 , 2 J2 , . . . , V h } . 

Theorem 5.4. Let q = p , for p an odd prime and t a non-negative integer, and let n = r ■ q, 
where r is a positive integer. Suppose that there is a nonessential element j G [n — 1] that is 
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not divisible by q. Furthermore, suppose that there exist integers a and b not divisible by p such 
that a = b mod 2, and (3 n {S) is congruent to either a or b modulo p for all SC [n — 1]. Then 

\{S C [n - 1] : (3 n (S) = a mod 2p}\ = \{S C [n - 1] : p n (S) = b mod 2p}\, 

\{SC [ n -l] : f3 n {S) = a + p mod 2p}\ = \{S C [n - 1] : f3 n (S) = b + p mod 2p}\, 

In the case when the proportion p(n) is 1/2, the four cardinalities above are all equal to 2 n ~ 3 . 

Proof. Consider the collection of sets S C [n — 1] such that (3 n {S) = a = b mod 2. For S 
in this collection such that j $ S, we have (3 n (S) = {3 n (S U {j}) mod 2. However, since g 
does not divide j, we have /3 n (S) = —0 n {S U {j}) modp by Proposition 15.1) that is, (3 n (S) ^ 
U {j}) modp, as a (resp., b) is not congruent to —a (resp., —6) modulo p. Hence this 
collection splits into two classes of equal size when divided according to the value of (3(S) 
modulo 2p. The same argument holds for the sets S satisfying j3(S) = a + p = b + p mod 2. □ 

By the Chinese remainder theorem, we have (3(S) = ±l,p ± 1 mod 2p. For non-Mcrsennc 
primes we can say more. (Recall that p(n) is defined as the ratio of the number of subsets 
SC. [n — 1] such that fi n {S) is odd to the total number 2™ _1 of subsets of [n — 1].) 

Theorem 5.5. Let q = p l be an odd prime power which has k 1 's in its binary expansion. 
Suppose that q > 2 k — 1, that is, q is not a Mersenne prime. Then 

\{SC [q-l] : (3 q (S) = I mod 2p}\ = \{S C [q - 1] : q (S) = -1 mod 2p}\, 

|{5C[j-l] : P q (S)=p-lmod2p}\ = |{SC[ g -l] : /3 q (S) = p + 1 mod 2p}\, 

In the case the proportion p(q) is 1/2, the four cardinalities above are equal to 2 q ~ 3 . 

Proof. Since q > 2 k — 1, as in the proof of Theorem 13.41 there exists a nonessential clement 
j E [q-l]- Thus Theorem IQ1 applies. □ 

When q = p is a prime and k — 2, Theorem 15.51 applies only to the Format primes which 
are greater than 3, that is, 5, 17, 257 and 65537. We also know that the proportion is 1/2 for 
the case k = 3, that is, primes whose binary expansion has three l's. The first few such primes 
are 7, 11, 13, 19, 37, 41, 67, 73, 97; see sequence A081091 in The On-Line Encyclopedia of Integer 
Sequences. 

For prime powers of the form q = p* with t > 2, the only case with k = 2 we know is 
q = 3 2 . Similarly with k = 3 we know six cases: 5 2 , 7 2 , 3 4 , 17 2 , 23 2 , 257 2 and 65537 2 . It is 
not surprising that the squares of the Fermat's primes and the square of 3 2 appear in this list. 
The two sporadic cases are 7 2 and 23 2 . 

The next theorem concerns permutations of size twice a prime power. 

Theorem 5.6. Let q = p be an odd prime power which has k 1 's in its binary expansion. 
Then 

\{S C [2q - 1] : !3 2q {S) = 1 mod 2p}\ = \{S C [2q - 1] : (3 2q (S) = -1 mod 2p}|, 
|{5C \2q-l] : /3 2(? (5)=p-lmod2p}| = \{S C [2q - 1] : (3 2q (S) = p + 1 mod 2p}|, 
In the case the proportion p(q) is 1/2, the four cardinalities above are equal to 2 2q ~ 3 . 
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Proof. By Corollary 15. 3( (32 q (S) = (— l)' 5 ' q }> mod p. Furthermore, since 2q is even, the ele- 
ment 1 is nonessential, and Theorem 15.41 applies. □ 



6. The descent set polynomial 

A different approach to view the results from the previous sections is in terms of the descent 
set polynomial 

Qn(t)= J2 t0n(S) - 
SC[n-l] 

The degree of this polynomial is the ?ith Euler number E n . For n > 2 the polynomial is divisible 
by 2t. Theorems 12 . 1 [ [5751 15.51 and 15.61 can be reformulated as follows. 

Theorem 6.1. (i) For a positive integer n we have Q n {— 1) = 2™ ■ (1/2 — p{n)). In partic- 
ular, when n has two or three 1 's in its binary expansion, then —I is a root of Q„ (t) . 

(ii) For n — 2 J > 4 the imaginary unit i is a root of Q n {t). 

(iii) Let q = p be a prime power, where p is an odd prime. Suppose that q has k 1 's in 
its binary expansion and satisfies q > 2 — 1. Let £ be a primitive 2p-th root of unity. 
Then 

Q q (0=2*-Re(0-(p(q)-± 

where Re(C) denotes the real part of Q. 

(iv) Let q — p be a prime power, where p is an odd prime. Suppose that q has k l's in its 
binary expansion. Let £ be a primitive 2p-th root of unity. Then 



) 2g (C)=2 2 «-Re(C)-(p(g)-0 



It is curious to observe that the polynomial Q n (t) quite often has zeroes occurring at roots 
of unity. An equivalent formulation is that Q n (t) often has cyclotomic polynomials $>k(t) as 
factors. (Recall that the cyclotomic polynomial &k(t) is defined as the product 0^(^ — C)i where 
£ ranges over all primitive fcth roots of unity.) See Table [2] for the cyclotomic factors of Q n {t) 
for n < 23. 

Lemma 6.2. Let q be an odd prime power. Then the cyclotomic polynomial <5> q does not divide 
the descent set polynomial Q n (t). 

Proof. If q is a power of an odd prime p, then <E> g (l) = p. Since Q n {l) = 2 n ~ 1 has no odd 
factors, the lemma follows. □ 

Lemma 6.3. Let q be the odd prime power p 1 . Then 

(i) If n = 2 3 then the cyclotomic polynomial $2g does not divide Q n {t). 

(ii) If n has four 1 's in its binary expansion and p > 5 then the cyclotomic polynomial $>2q 
does not divide Q n (t). 

(iii) If n has five 1 's in its binary expansion and p > 11 then the cyclotomic polynomial <I>2 g 
does not divide Q n (t). 
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Proof. We have 3>2q(— 1) = P- Since Q n {— 1) = 2™ • (1/2 — p(n)) the result follows by consulting 
Table CD □ 



7. Quadratic factors in the descent set polynomial 

In order to study the double root behavior of the descent set polynomial Q„ (t) or, equiva- 
lently, quadratic factors in Q n (t), we need to prove a few identities for the descent set statistics. 
We begin by introducing the multivariate ab- and cd-indexes. Let ai, a2, . . . and bi,b2, . . . be 
non-commutative variables. For S C [n — 1] let us be the monomial u\U2 ■ ■ ■ !i n _i where Ui = &i 
if i G" S and Ui = hi if i G S. The multivariate ab-index of a poset P of rank n is defined as 
the sum 

V{P) = J2 h s-us, 

s 

where S ranges over all subsets of [n— 1]. 

Lemma 7.1. For an Eulerian poset P the multivariate ab-index can be written in terms of the 
non-commutative variables Ci = a.; + b^ and dij + ± = ajb^ + i + bjaj+i. 



Proof. Observe that by adding the index i to the ith letter in an ab-monomial of degree n— 1, 
we obtain a natural bijection between the regular and the multivariate ab-indices of the same 
poset P. Thus the statement of the lemma is equivalent to the statement that the regular 
ab-index of P can be expressed in terms of the variables c = a + b and d = ab + ba. □ 



In this case, we call the resulting polynomial the multivariate cd-index. Observe that for 
a rank n Eulerian poset each of the indices 1 through n appears in each monomial of the 
multivariate cd-index. 

Proposition 7.2. Let hs be the flag h-vector of an Eulerian poset P of rank n, or more 
generally, hg belongs to the generalized Dehn-Sommerville subspace. Let T C [n — 1] such that 
T contains an interval [s, t] = {s, s + 1, . . . , t} of odd cardinality with s — l,t + 1 £ T . Then 

J2 (-i) |snT| -/ ls = o. 

SC[«-1] 



Proof. The sum is obtained from the multivariate ab-index of the poset P by setting a.; = 1 
and 

-l if i eT, 

1 otherwise. 

Notice that Ci = for i G [s, t] and that d s _i iS = d tt t+i =0. If s = 1 we set do,i — 0, and if 
t = n — 1 we set d n _i .„ = 0. Since P is Eulerian, the multivariate ab-index can be written in 
terms of multivariate cd- monomials (Lemma 1 7. 1|) . A multivariate cd- monomial that contains 
d s _i iS or d tit+ i evaluates to zero. Since the interval [s,t] has odd size, a multivariate cd- 
monomial not containing d s _i iS and d tjt+ i must contain at least one variable with i G [s, t\. 
Hence this monomial also evaluates to zero. □ 
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Observe that the identity in Proposition 17.21 is a part of the generalized Dehn-Sommerville 
relations; see pQ. 

Theorem 7.3. If the binary expansion of n has two 1 's and n > 3, then "J" 2 , divides Q n (t). 



Proof. Suppose that n = mi + m 2 , where m\ = 2 n , m 2 = 2 32 , and ji > j 2 . From the proof of 
Theorem 13.41 we have 

1 mod 2 if \Sn{mi,m 2 }\ = 0,2, 
mod 2 if \Sn{m 1 ,m 2 }\ = 1. 



/?„(S) 
Hence 



q;(-i) = E &(*)-(-i) /s " {s) - 1 

SC[ra-l] 

E (-l) |sn{mi,ma}l -/3„(5), 

SC[n-l] 

which is zero by Proposition [721 D 
Theorem 7.4. If n = 2 j > 4 i/ien $4 divides Q n {t). 



Proof. Let to = n/2. The proof of Theorem EH states that f3(S) = (-l)\ s -t m }\ mod 4. Let i 
be the imaginary unit, so that i 2 = —1. Observe that i^ -1 ) _1 = (— l) fc . We have 

QUi) = £ f3 n (S).i^ s ^ 

SC[n-l] 

= E (-l) |5 - {m}l -/3«(S). 

SC[n-l] 

By Proposition [7J21 QnW = 0, since S — {to} = «S fl {1, . . . , to — 1, to + 1, . . . , n — 1} and to — 1 
is odd. □ 



The next result applies to prime powers that have two l's in their binary expansion. The 
only cases known so far are the five known Fermat primes 3, 5, 17, 257, 65537 and the prime 
power 3 2 . 

Theorem 7.5. Let q = p l be a prime power, where p is an odd prime and assume that q has 
two l's in its binary expansion. Then the cyclotomic polynomial divides Q 2q {t). 



Proof. In this case n = 2q = m + 2, where m = 2 3 . From the proof of Theorem 13 .41 we have 

, ,f 1 mod 2 if |Sn {2,m}| = 0,2, 
P 2q (b)-< Q mod2 if |5- n {2,77x}| = 1. 



Hence combining it with the proof of Corollarv l5.31 we have 

(_l)|S-{g}| mo d2p if \SC\ {2,to}| = 0,2, 



MS) = 



P 



+ (_l)|S-{g}| mo d2p if |S*n {2,to}| = 1. 
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Thus for £ = Re(£) + Im(£) • i a 2p-th primitive root of unity, we have that 

= (_l)|Sn{2,rn}| . ^ Re(c) + . Im(C ) . A . 

Evaluating the sum and using the fact that 

\S n {2, m}\ + \S - {q}\ = \Sf) {2, q, m}\ mod 2 

we have 

C-Q 29 (C) = E W5)-C W5) 

SC[2g-l] 

= Re(C)- 2 (-l) |Sn{2 ' m}L /3 29 (5) 

SfC[2Q[-l] 

+ Im(C)-i- £ (-l)' Sn ^™>'-/3 29 (5), 

SC[2g-l] 

where both sums vanish by Proposition [721 D 

8. The signed descent set polynomial 
Similarly to the descent set polynomial we can define the signed descent set polynomial: 

Qt(t) = E ^ ±(s) ; 

SC[n] 

see Section U for definitions relevant to signed permutations. The degree of this polynomial is 
the nth signed Euler number , which is the number of alternating signed permutations of 
size n. Yet again, for n > 1 this polynomial is divisible by 2t. Theorem 14.61 can now be stated 
as follows. 

Theorem 8.1. Forn> 2 the signed descent set polynomial (t) has the cyclotomic factor $4. 

The space of quasisymmctric functions of type B is defined as BQSym = Z[s] QSym. 
Quasisymmetric functions of type B were first defined by Chow [5]. We will view them to be 
functions in the variables s, Wi,W2, ■ ■ ■, that are quasisymmetric in u>i, u> 2 , . . .. For a composition 
(70, . . . , 7 m ) define the monomial quasisymmetric function of type B by 

M (70,T 1 ,..,T m ) = s7 °" 1 - M (7 1 ,...,T m )- 

A third method to encode the flag vector data of a poset P of rank at least 1 is the quasi- 
symmetric function of type B 

(8-1) F B {P)=Y,Mf {cy 

c 

where the sum is over all chains c = {0 = xq < xi < ■ ■ ■ < x m = 1} in the poset P; see [12]. A 
different way to write equation (|8.ip is 

(8.2) F B (P)= E s"^- 1 ■ F([x,l}). 

6<x<i 
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The diamond product of two posets P and Q is 

PoQ = (p-{6})x(Q-{6})u{6}. 



Using identity (|8.2|) one can show that the type B quasisymmetric function of a poset is multi- 
plicative with respect to the diamond product of posets, that is, Fb(PoQ) = Fb(P) ■ Fb(Q)- 
Applying the bijection D between compositions and subsets, we have 

f b(P)= E fs-Mi, 

SC[n-l] 

where we write Mg instead of M^, s ^, and the poset P has rank n. The fundamental quasi- 
symmetric function of type B , denoted by L^, is given by 

14 = y, Ml 

TCS 



Then the flag h-vector appears as the coefficients in the decomposition 

F B (P)= J2 h S' L l 

SC[n-l] 

where the poset P has rank n. 

The cubical lattice C n has rank n + 1 and can obtained as a diamond power of the Boolean 
algebra B2, that is, C n = B^ 1 . Therefore we have the following result. 

Lemma 8.2. The type B quasisymmetric function of the cubical lattice is given by 

F B (C n ) = (s + 2-M {1) ) n . 

Theorem 8.3. For p an odd prime the cyclotomic polynomial <&4 P divides the signed descent 
set polynomial Qp(t). 

Proof. Observe that modulo 4 we have 

F B (C P ) = (s + 2-M (1) )P 

= sP + 2-p-sP- 1 -M (1) 

= M ( P +i) + 2 ■ P ■ M ( P ,i) 
= M + 2 • M {p} 

= Y(-l) lSl -L B s +2-Y(-l) lSl - 1 -L B s 

sc\p] P es 

- EM^if + EM^-Lf mod4. 

p<ZS pes 

Hence the signed descent set statistics satisfy /3p(S) = (— 1)I s ' _ {p} mod 4 for S C [p\. Now 
modulo p we have 

F B (C P ) = (s + 2-M (1) ) p 
= S P + 2-M {p) 
ee M* +x) +2.Ml p) 
ee Mg + 2 ■ Mm mod p. 
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(8.3) 



This directly implies that Pp(S) = (— l)' s mod p. Combining these two statements wc 
obtain 

(-l) |s| mod4p ifl,p£S, 

(-l)^- 1 mod4p ifl,pe5, 

2-p+(-l)l s ' mod4p HlgS,peS, 

2 -p+ (-l)! 5 !" 1 mod4p ifl€S,p#S. 

Observe that for p <£ S we have 0p(S) = fif(S U {p}) + 2 • p mod 4p, implying that ^ = 

C ,Su W) for C a 4p-th primitive root of unity. Now sum over all subsets of [p], and the result 
follows. □ 

Theorem 8.4. For p an odd prime, $| does not divide the signed descent set polynomial 
Qp (t) . In fact, evaluating the derivative of the signed descent set polynomial (t) at (, where 
£ is a Ap-th primitive root of unity, gives 

C ■ Qf(C) = Im(C) • i • (-l)^ 1 )/ 2 -V-p- E p -t. 



S-{1} 



Proof. From (|8. 3[) we have: 

SC[p] 

= £ /? P ± (^)-(-i) |sn{1 ^ }l -c ( - 1 

SC[p] 

= Rc(C)- £/£($)• (-l)^* 1 ^ 

+i m (o-i-E^)-(- 1 ) |sn{1,p}L (- 1 ) l ^ {1}l - 

The first sum is zero by Proposition [721 The second sum simplifies to 

MO-i- E/^SM-i) 1 ^- 111 . 

SQIp] 

This sum can evaluated by setting a.j = 1, bi = • ■ ■ = b p _i = — 
multivariate ab- index of the cubical lattice C p . Observe that Ci 



1 and bp = 1 in the the 



Cp-l 



and 



dp-i, p = 0. Hence the only surviving cd-monomial is dx,2 • ■ ■ d p _2. p -ic p . The coefficient of 
this monomial is computed as follows: 

d (p-i)/2 c l 9 ( C j = 2 (f" 1 )/ 2 . [(2d)( p - 1 )/ 2 cl *(C p ) 



2 (p-l)/2 



(ab) 



(P"l)/2, 



a-*(B p 



( ab )(p-i)/2 b l a . ^( Bp ) 

= 2(f- 1 )/ 2 -p- [b(ab)^- 3 )/ 2 " 
= 2^^-p-E^. 



*(^ P -i) 



The third step is MacMahon's "Multiplication Theorem"; see [131 Article 159]. It can be stated 
in terms of the ab-indiccs as follows: 

' m + 
m 



[uav]^(B m+n ) + [ubv]m{B. 



rn+n J 



[u]^(B m ) + [v]*(B n ), 
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where u and v have degrees m — 1 and n — 1, respectively. The monomial itself evaluates to 
(_2)(p~ 1 )/ 2 . 2, since dj.,2 = • • ■ = d p _2,p-i = —2 and c p = 2. Combining all the factors, the 
evaluation at £ follows. □ 



n 


degree 


cyclotomic factors of Q n 


(t) 


3 


2 


<j> 2 








4 


5 


*! 








5 


16 


*! 








6 


61 


*i 


*i • $io 






7 

8 


272 
1385 


<j> 2 

*! 


$28 






9 


7936 


*i 


*6 • $18 






10 


50521 


*l 


$6 • 


• $18 • $30 




11 


353792 




$6 ' *22 






12 


2702765 


*i 


$6 ' $10 


• $18 ' $22 ' 


$22 ' $66 ' $110 ' $198 


13 


22368256 




*26 






14 


1.993- 10 8 


*2 


$ 2 ' $4 ' 


*14 ' $14 ' 


$26 ' $28 ' $182 


15 


1.904- 10 9 










16 


1.939 • 10 10 


*2 


$12 ' $2f 


• $44 ' $52 


' $60 ' $156 ' $220 ' $260 


17 


2.099- 10 11 


*l 








18 


2.405 • 10 12 


*l 


*i • $i8 


• $34 ' $102 


• $306 


19 


2.909 • 10 13 




*38 






20 


3.704 • 10 14 


*i 


$6 ' $10 


• $30 ' $34 ' 


$1 8 • $102 ' $114 ' $170 






•$190 • $510 ■ 


3*570 ' $646 


$1938 ' $3230 ' $9690 


21 


4.951 • 10 15 


$2 


$6 ' $14 


■ $42 




22 


6.935- 10 16 


$2 


$2 ' $14 


• *22 ' $22 


$154 


23 


1.015- 10 18 











$5 



Table 2. Cyclotomic factors of Q n (t). 



9. Concluding remarks 



Is there a reason why p(n) — 1/2 factors so nicely? See Tabic Q] 

The two main results for unsigned permutations in Section [21 Theorems 12.11 and 13 . 4[ can also 
be proved using the ab-index and the mixing operator; see [7j. We have omitted this approach 
since Kummcr's theorem and the quasisymmetric functions are more succinct in this case. 

Tables [5] and [3] contain cyclotomic factors of polynomials Q n {t) and Qn{t) for small n. Those 
factors whose presence is explained in this paper are highlighted in boldface. Here are several 
observations about the data in Tabled 

(i) All the indices k of cyclotomic factors $& of the polynomials Q n (t) are even. 

(ii) Any prime factor p that occurs in an index of a cyclotomic factor of Q n (t) is less than 
or equal to n. 
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n 


degree 


cyclotomic factors of Q^(t) 




2 


3 


$4 








3 


11 


$4 


■ $8 ' 


$12 




4 


57 


$4 


■ $16 


• $32 




5 


361 


$4 


■ $16 


• $20 ' $32 ' $80 




6 


2763 


$4 


. $ g . 


$24 ' $32 • $40 • $96 • 


$120 • 


7 


24611 


$4 


■ $8 ■ 


$24 • $28 ' $32 ' $56 


■ $168 ' 


8 


250737 


$4 


■ $32 


' $64 ■ $224 • $448 ' $ 


512 


9 


2873041 


*4 


. $ 12 


• $32 • $36 ■ $64 ■ $96 


• $192 



'$512 ' $576 ' $1344 ■ $1536 ■ $4032 1 $4608 

10 36581523 * 4 ■ $8 ■ $24 ■ $40 • $64 ■ $72 • $120 • $192 ■ $320 ■ $360 

•$448 ' $512 • $960 ' $1344 ' $1536 " $2240 " $2560 " $6720 ' $7680 

11 5.123 • 10 8 * 4 ' $8 ' $40 ' $44 ' $64 ' $88 ' $192 ' $320 ' $440 

•$512 • $704 • $960 ' $2112 • $2560 - $3520 • $5632 

12 7.828 • 10 9 $4 ' $16 ' $32 ■ $48 ' $96 " $160 " $176 ■ $288 ' $352 

•$480 ' $512 • $528 ' $1056 ' $1440 ' $1536 ' $1760 • $2560 
•$3168 ' $4608 ' $5280 - $5632 

13 1.296 • 10 11 * 4 ' $16 ' $32 ■ $48 ' $52 ' $160 ' $208 ' $352 " $416 

•$512 ' $624 ' $ 1536 ' $1760 - $2080 ' $4576 ' $5632 ' $6656 

14 2.310 • 10 12 $ 4 ' $8 ■ $32 • $56 • $104 " $224 " $352 ■ $416 ' $512 

•$728 • $1536 - $2464 ' $2912 ' $3584 • $4576 " $5632 ' $6656 

15 4.411 • 10 13 * 4 ' $8 ' $12 ' $20 • $24 ' $32 • $40 ' $56 ' $60 

•$96 • $120 ' $160 ' $168 - $224 ' $280 ' $416 ' $480 ' $512 

•$672 ' $840 ' $1120 ' $1248 " $1536 ' $2080 ' $2560 ' $2912 • $3360 

•$5632 ' $6240 ' $6656 • $7680 - $8736 

16 8.986 • 10 14 * 4 ' $64 ' $ 128 • $192 • $320 ' $640 ' $896 ' $960 

•$1024 • $1664 • $3072 • $4480 • $5120 ' $8320 

17 1.945 • 10 16 * 4 ' $64 ' $68 ' $128 ' $640 ' $896 

'$1024 ' $1088 ' $2176 1 $4480 " $5120 

18 4.458 • 10 17 * 4 • $8 ■ $24 • $72 • $128 " $136 " $384 ■ $408 ' $640 

•$1024 ' $1152 ' $1224 ' $1920 " $2176 ' $3072 1 $5760 - $6528 - $9216 

Table 3. Cyclotomic factors of Qt(t)- 



(iii) If $fe x and $fe 2 arc factors of Q n (i), so is $ gc d(fci,fc 2 )- That is, the set of indices is closed 
under the meet operation in the divisor lattice. 

(iv) If hi divides k2, &2 divides and $^ 1 and $fc 3 occur as factors in Q n (t), then so does 
$fc 2 . This is convexity in the divisor lattice. 

(v) If both $fcj and $fc 2 divide Q n (t), where fei divides k2, then the multiplicity of $fc x is 
greater than or equal to the multiplicity of $£ 2 . 

(vi) If p is not a Mersenne prime then the largest cyclotomic factor occurring in Q p (t) is 

$2p- 

(vii) When p(n) ^ 1/2 then there are no cyclotomic factors in the descent set polynomial 
Qn(t). 

(viii) For all primes p we conjecture that $ 2 p divides Q2p- 



Moreover for the signed descent set polynomial we observe that: 
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(ix) For n > 3 the the cyclotomic polynomial <E>4„ divides the signed descent set polynomial 

Qt(t)- 

(x) For n > 5 the the cyclotomic polynomial *&4 n t n -i) divides the signed descent set poly- 
nomial Qn(t)- 

Can these phenomena be explained? 

For what pairs of an integer n and a prime number p does the descent set statistic (3 n {S) 
only take two values modulo pi 

Finally wc end with two number-theoretic questions. Are there infinitely many primes whose 
binary expansion has three l's? The only reference for these primes we found is The On-Line 
Encyclopedia of Integer Sequences, sequence A081091. Are there any more prime powers with 
two or three ones in its binary expansion? 
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